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We consider a novel class of hard branching hadronic processes a + b — > c + d + e, where hadrons c 
and d have large and nearly opposite transverse momenta and large invariant energy which is a finite 
fraction of the total invariant energy. We use color transparency logic to argue that these processes 
can be used to study quark generalized parton distributions (GPDs) for baryons and mesons in 
hadron collisions hence complementing and adding to the studies of GPDs in the exclusive DIS 
(deep inelastic scattering) processes. We propose that a number of GPDs can be investigated in 
hadron facilities such as J-PARC (Japan Proton Accelerator Research Complex) facility and GSI- 
FAIR (Gesellschaft fur Schwerionenforschung -Facility for Antiproton and Ion Research) project. 
In this work, the GPDs for the nucleon and for the N — > A transition are studied in the reaction 
N + N —> N + 7r + B, where N, tt, and B are a nucleon, a pion, and a baryon (nucleon or 
A), respectively, with a large momentum transfer between B (or n) and the incident nucleon. In 
particular, the Efremov-Radyushkin-Brodsky-Lepage (ERBL) region of the GPDs can be measured 
in such exclusive reactions. We estimate the cross section of the processes 7V + iV^iV + 7r + i?by 
using current models for relevant GPDs and information about large angle ttN reactions. We find 
that it will be feasible to measure these cross sections at the high-energy hadron facilities, and get 
novel information about the nucleon structure, for example, contributions of quark orbital angular 
momenta to the nucleon spin. The studies of N —> A transition GPDs could be valuable also for 
investigating electromagnetic properties of the transition. 

PACS numbers: 13.85.-t, 13.60.Lc, 12.38.-t 



I. INTRODUCTION 

Understanding complexity of the structure of hadrons 
is one of prime objectives of QCD (Quantum Chromo- 
dynamics). For a long time, the main focus was the 
study of single parton densities in inclusive processes for 
which QCD factorization theorem is valid More re- 
cently, many investigations started to focus on nucleon 
spin structure by polarized lepton-nucleon scattering and 
proton-proton collisions. It became clear that only a 
small fraction (20—30 %) of nucleon spin is carried by 
quarks 0] . There are still a large uncertainty in the gluon 
spin contribution however, recent measurements in- 
dicate a tendency that it is also small in the region of 
momentum fraction i ~ 0.1 j|. It is, therefore, impor- 
tant to find a contribution from orbital angular momenta 
as a possible source for explaining the nucleon spin. 

A parallel development of the last two decades was re- 
alization that the QCD factorization theorem is valid for 
deeply virtual Compton scattering Q and for exclusive 
production of mesons in deep inelastic scattering (DIS) 
by longitudinally polarized photons 0, [f| . Amplitudes of 
these processes are expressed through universal objects 



- generalized parton distributions (GPDs) which contain 
information about both longitudinal and transverse dis- 
tributions of hadrons. 

In the forward scattering limit, several GPDs be- 
come usual parton distribution functions with the trans- 
verse momentum dependence connected to the trans- 
verse coordinate distribution of partons in the hadron [?J . 
Therefore, studies of GPDs allow to establish the three- 
dimensional image of the nucleon and hence provide crit- 
ical information for description of nucleon-nucleon colli- 
sions at collider energies [fl . The first moments of GPDs 
are given by the corresponding hadron form factors pro- 
viding a link between high-energy and medium-energy 
dynamics. At the same time the second moment of a 
certain combination of GPDs is related to the total an- 
gular momentum contribution to the nucleon spin and 
hence probes orbital angular momentum effects in the 
spin structure of nucleon Q . 

For the recent reviews of the properties of GPDs and 
applications to the study of the hard exclusive DIS pro- 
cesses, we refer the reader to Refs. [ToL [Tlj . Several 
theoretical calculations have been performed [HI, [H, 
EH fl5|. [l6L and we will use the ones described in 
Ref. [T2l • Experimental studies of GPDs are performed 



at HERA (Hadron-Elektron-Ringanlage) collider and in 
several fixed target experiments [181 Il9j | . It was also pro- 
posed to study GPDs in the irN — > + N process 

The QCD factorization theorem for exclusive hard pro- 
cesses induced by longitudinally polarized photons is 
valid for any exclusive final states in particular exclu- 
sive meson production [f|. Hence, one can probe various 
non-diagonal transitions between baryons - transitional 
GPDs III. For example, the GPDs for the N -> A tran- 
sition, where N and A indicate the nucleon and A(1232), 
should be important for the studies of nucleon spin [2lj ]. 
There is also an important link to the ./V — > A transi- 
tional electromagnetic form factors, which are studies in 
the low-energy electron experiments [22j|. Electromag- 
netic moments of unstable particles such as the A can- 
not be measured by usual methods for stable nuclei [23| . 
Electroma gnet ic moments for the N — ► A are interesting 
quantities [22|]. In particular, the C2 and E2 transition 
moments indicate shape deformations of the nucleon and 
A, so that they are valuable for pro ving interesting in- 
formation on quark-gluon dynamics 24(. The N — > A 
GPDs are related to the nucleonic GPDs according to a 
chiral-quark soliton model in the large N c limit [l(j, Hl[ , 
where N c is the number of colors. Using the relation to 
the nucleonic GPDs, we could learn about contributions 
of quark orbital angular momenta to the nucleon spin. 

It was suggested in Refs. (2|| [26| that one can in- 
vestigate presence of small-size color singlet qq and qqq 
clusters in hadrons using large-angle branching hadronic 
processes a + b — ► c + d + e where the hadron e is pro- 
duced in the fragmentation of b with fixed Feynman xf 
and fixed transverse momentum p^' , while the hadrons c 
and d are produced with large and near balancing trans- 
verse momenta: p T m —p T . 

In this paper, we will argue that investigations of these 
processes allow to probe various diagonal and transitional 
GPDs of hadrons. We will focus on one of the simplest 
subclass of these reactions N + N — ► N + it + B, where 
7r indicates the pion and B does the nucleon or A. We 
will argue that these reactions allow to probe the GPDs 
for the nucleon and the N — ► A transition in the kine- 
matical region of Efremov-Radyushkin-Brodsky-Lepage 
(ERBL) [27| at hadron facilities in addition to the lep- 
ton scattering measurements. Especially, there are fu- 
ture projects such as the J-PARC (Japan Proton Ac- 
celerator Research Complex) facility [28] and GSI-FAIR 
(Gesellschaft fur Schwerionenforschung -Facility for An- 
tiproton and Ion Research) project [2!| for investigat- 
ing high-energy proton and anti- proton reactions. Other 
branching reactions which one can study at these fa- 
cilities include reactions initiated by meson projectiles 
which provide a unique way to probe meson GPDs: 
7r + N — > 7T + 7T + low-pT baryon and processes initi- 
ated by protons (antiprotons) involving production of 
two large angle baryons (baryon and antibaryon at the 
GSI-FAIR) with a small transverse momentum for the 
pion: p(p) + N — *• p(p) + N + low-px meson. The cross 



sections of the processes with two large transverse mo- 
mentum baryons are expressed by skewed (transitional) 
parton distributions N — > M [26[, where M indicates a 
meson. Experimental indications for existence of back 
to back production of two protons with large pt were re- 
ported in Rcf. [30| . 

This article is organized as follows. In Sec. [Til 
we provide arguments for validity of factorization for 
two — > three process and introduce the description of 
the N + N — > N + tt + B cross section in terms of a 
N — > B amplitude and a meson- nucleon scattering part. 
The GPDs for the nucleon and the N — ► A transition 
are defined and their major properties are introduced in 
Sec. IIIIl The N — > B amplitude is expressed through 
the nucleonic or N — > A GPDs in Sec. IIVI The meson- 
nucleon scattering cross sections are parametrized so as 
to explain experimental measurements in Sec. [Vj The 
N + N^N + tt + B cross sections are shown in Sec. IVII 
by assuming a J-PARC kinematics for the initial proton 
and by taking a typical model for the GPDs. The results 
are summarized in Sec. IVIII 



II. FACTORIZATION OF CROSS SECTION FOR 

N + N^N + tt + B 

l I 

It was argued in Ref. [31| that cross sections of large 
angle scattering processes: a + b — > c + d in the limit of 
s — > oo with # cm =const, where s is the center-of-mass 
(cm.) energy squared and 9 cm is the scattering angle, 
are dominated by the short-distance configurations in 
hadrons both in the initial and final states and satisfy di- 
mensional counting rules. The data on a variety of such 
processes (see Ref. HH and references therein) are con- 
sistent with the dimensional counting rules. They also 
indicate that the processes, where quark exchange is al- 
lowed, are much larger than those where it is forbidden. 
The especially interesting is the observed systematics of 
reactions induced by mesons: the ratio of K + p — > K + p 
and it + p — > tt + p cross sections at 9 cm = 90°, where p is 
the proton, is close to the square of the wave functions 
K + and tt + mesons at the origin: [/if/j^] 2 ~ 1-4, while 
it is ~ 0.5 for a small momentum transfer \t\. 

A test of dominance of the small-size configurations in 
these processes was suggested by Brodsky and Mueller 
[H,[3J]. They pointed out that since small-size configura- 
tions weakly interact with other hadrons the cross section 
of the large angle processes like a + A—t a + p+ (A — 1) , 
where a and A indicate a hadron and a nucleus with mass 
number A, should be proportional to the number of pro- 
tons in the nucleus. This regime is usually referred to 
as the color transparency (CT). The phenomenon which 
slows down the onset of the CT regime is the space-time 
evolution of the small wave packages. Evolution occurs 
over the distance scale l co h (coherence length) , which can 
be estimated as l co h ~ 0.3 -j- 0.4 fm • ph/GeV, where ph 
is the hadron momentum [35| . Recently, the color trans- 
parency phenomenon was observed in the electroproduc- 
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tion of pions [36| with the change of the transparency 
consistent with the prediction of Ref . [13] , in which above 
estimate of l co h was used. For the review of the CT phe- 
nomena in the high-energy processes, see Ref. Q . 

Overall, it appears that there is a good evidence for 
dominance of the small-size configurations in the large 
angle mcson-nucleon scattering and for suppression of 
the interaction of small-size configurations near the in- 
teraction point. Hence, the mcson-nucleon scattering at 
large angles should be a good probe of short-distance 
effects already at moderate energies s-^n ~ 12 GeV 2 
and large 9 cm ~ 90°. Therefore, we start our investi- 
gation of the two — ► three processes with the reaction 
N + N^N + ir + B in the kinematics where nucleoli 
scatters at large \t'\ off a small-size color singlet qq clus- 
ters in the nuclcon. This is analogous to the process of 
7* + TV — > 7r + B scattering at x > 0.2, where dominant 
contribution originates from the scattering off small-size 
color singlet qq clusters in the nucleon [H|. These pro- 
cesses allow to investigate the nuclconic and N — > A 
transition GPDs. An additional advantage is that these 
GPDs are extensively studied in the literature so that we 
will be able to use existing analyses of these GPDs to 
perform numerical estimates. In this section, we explain 
how to calculate the cross section of these reactions and 
give arguments for validity of the factorization. 

The cross section is given by 



da 



S 



1 d 3 p c 1 d 3 p d 1 d 3 p e 
X 2E C (2tt) 3 2E d (2tt) 3 2E e (2tt) 3 

X (27r) 4 5 4 (p a + p b - p c -Pd~ Pe), 



NN—*NttB 



(1) 



where Mnn->NttB is the matrix element for the process 
NN — > NnB, mjv is the nucleon mass, and A^ (i=a, b, d, 
e) is a spin of an initial or final particle. The factor S is 
1/2 if the final baryon is identical to the final nucleon N 
and it is 1 if they arc not identical. The particle indices a, 
&, c, d, and e are assigned to the initial nucleons and final 
nucleon, pion, and baryon B (nucleon or A) as shown in 
Fig. [T] Their momenta and energies are denoted as pi 
and Ei (i=a, b, c, d, e). Averages arc taken for the initial 



N (d) 




spins A Q and A&. 
defined 



Mandelstam variables s. t, s', and t' are 



S = {Pa +Pb) , t =(p a - p e ) , 
s' = (Pc + Pd) 2 , t'=(p b -p d ) 2 , 



(2) 



for this reaction. Note that two complementary kine- 
matic regions contribute to the discussed cross section. 
In one kinematics B has small momentum relative to nu- 
cleon (a) - the target nucleon, while in the second kine- 
matics t is small relative to the projectile - nucleon (b). 
These two cross sections are obviously equal. The cor- 
responding amplitudes do not interfere in the discussed 
limit. Hence for simplicity wc will present the results for 
one kinematics only. At the moment, there exists no pro- 
cedure for calculating cross sections of N+N — > N+n+B 
process in a generic kinematics. However, it is simplified 
for the hard regime: \t'\ ^> m 2 N and \u'\ = \{pb — p c ) 2 \ ^> 



m 



N 



Formal arguments can be developed for the limit 



oo, — = const, 
s' 



(3) 



since the leading-order QCD diagrams dominate the cross 
section of two-body processes [31( • All the partons which 
experience a large change of the direction have to be at- 
tached to hard lines. This idea is used as a counting 
rule for estimating energy dependence of exclusive cross 
sections at high energies [3l[|38j. In our studies, the vari- 
able \t\ is fixed and sufficiently small (\t\ <C mjy), and the 
center-of-mass energy is >/s ~ \/6Q VTOO GeV at the 
J-PARC. 

Similar to the logic of the QCD counting rules of Ref. 
[3l| , we separate diagrams with minimal number of hard 
gluons. The hard part of these diagrams has the same 
topologic structure as the diagrams for the corresponding 
elementary reaction M+N — » M'+B. A typical diagram 
of this kind is presented in Fig. [2l 

It is easy to check that other processes which con- 
tribute to the N + N — * N + n + B cross section, are 




FIG. 1: iY + iV^TV + TT + B reaction. 



FIG. 2: (Color online) A typical leading subprocess contribu- 
tion, which corresponds to Fig. [4] to the N + N — > N + n + B 
cross section. The spirals indicate gluon exchanges. The thick 
spirals show hard gluon exchanges, whereas the narrow one 
shows a soft gluon exchange. In the large \t'\ and \u'\ process, 
the quarks in the final baryon and pion should be attached to 
the hard gluon lines. 
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FIG. 3: (Color online) A typical subleading contribution to 
the N + N^N + n + A cross section. As explained in the 
text, the nucleon (a) splits as N — > qq + B* — > n + B, then 
the qq pair interacts with another nucleon (b) (N + qq — > AT). 
Notations for soft and hard gluons are the same in Fig. [2] 

suppressed by powers of s' in the limit of Eq. ([3J. For 
example, in the case of a diagram shown in Fig. [3J the 
process proceeds via N ^ qq + B* qq + n + B, where 
B* indicates an intermediate baryonic state, in the lower 
blob of this figure and subsequent transition TV + qq — > N 
in the upper one. In this case, the amplitude is small be- 
cause the qq propagation is suppressed by the factor l/\t'\ 
(\t'\ S> wtjy), whereas the corresponding factor in Fig. [2] 
is l/\t\ (\t\ -C m 2 N ). Furthermore, the large momentum 
(\t'\) is also involved in the lower blob (N — > qq + B*) in 
addition to the upper one. It means that more hard gluon 
exchanges are involved in this case, which suppresses the 
contribution in Fig. [3J 

Next, we need to demonstrate that additional soft in- 
teractions, which are not suppressed by powers of t', are 
canceled out. The easiest way is to consider the pro- 
cess in the rest frame of fragmenting nucleon. Neglecting 
space-time evolution of the wave packets, which is legit- 
imate in the limit of Eq. ^ , we observe that transverse 
size of the projectile and two outgoing hadrons near the 
interaction point is given by ~ l/y[FJ. Hence, the color 
transparency arguments essential for the proof of factor- 
ization for exclusive hadron production in DIS Q are 
applicable, and the soft interactions are suppressed by at 
least a power of 1/t'. Hence, we find that the leading 
contribution to the N + N^N + tt + B amplitude is 
given by the factorized form. 

The transition amplitude Mnn^NttB is written as a 
convolution of generalized parton distribution G(N — > 
B) of qq pairs in the nucleon with a hard interaction blob 
H of qq - Sq elastic scattering and qq (3q) wave functions 
of the initial nucleon (b), the final nucleon (d), and (c): 

M G(N -»• B) g> O ® ip f N . (4) 

Here, (£> denotes convolution of the different blocks, 
namely an integration over the light-cone fractions of 
constituents involved in the process, etc. There is also 
an implicit summation over the intermediate states. We 
can also rewrite Eq. ([4]) in somewhat more transparent 
form as: 

M-NN~>NirB = ^N~*hB ^hN^irN, (5) 



(a) 




FIG. 4: Factorization ofN + N^N + ir + B process into 
N — > h + B and h + N -> nN amplitudes. 

where h indicates the intermediate qq hadronic state and 
the nucleon, respectively. This corresponds to the dia- 
gram of Fig. 2J 

Consequently, the cross section has the following 
generic factorized form in the limit of Eq. ([3]): 

~, , 2 ^ JO = f( a 'PBT)4>{s\ 9cm), (6) 

where a and psr arc the light-cone fraction and trans- 
verse momentum carried by B (in the GPD notation 
a = (1 — 0/(1 + 0: see section 3). The cross section is 
separated into the function / (and cj)) which depends on 
the variables a and pbt (s' and cm ). The invariant en- 
ergy of the subprocess is s' = (1 — a)s, and the function <p 
is the cross section for the qq-3q scattering, and it is given 
by the quark counting rules as <p(s',O cm ) ~ (s') 7(^ CTO ), 
where "f(9 C m) is a function which depends only on the an- 
gle 9 cm , and n is the total number of all interacting ele- 
mentary fields. Phenomcnologically, it appears more nat- 
ural to use for numerical studies of experimental s', t' de- 
pendence of the elementary reaction ttN — > MN , where 
M is a meson, rather than the asymptotic quark counting 
values. 

Note the cross section as given by Eq. ^ does not 
depend on the direction of the transverse component of 
the momentum of B, or equivalently the cross section 
does not depend on the angle between the plane formed 
by B and the target nucleon and the plane of the hard 
reaction. This provides an important test of factorization 
as the rescatterings, which are discussed below, lead to 
break down of this factorization property. 

To determine at what t' we may expect the onset of 
factorization, we notice that the minimal condition is 
that qqq and qq wave packets, which are small in the 
interaction point, involved in the hard interaction block 
should remain small while passing by the spectator sys- 
tem. Hence, the distance over which a small wave pack- 
age expands or contracts - l co h, should satisfy the condi- 
tion 

hoh > r N - 0.8 fm, (7) 

corresponding to expansion of the outgoing hadrons and 
collapse of incoming proton outside the target proton 
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FIG. 5: Time evolution of wave packets of interacting hadrons 
in the 2^3 process. 



as shown in Fig. [5] Here, 0.8 fm indicates the size of 
the nucleon. A similar condition exists in the discussion 
of the color transparency phenomena with nuclei where 
it is necessary to have l co h > 2Ra, where Ra is a nu- 
clear radius, to reach the regime of complete color trans- 
parency. The studies of the evolution rates of qqq and qq 
wave packets in the framework of the nuclear color trans- 
parency have led to the estimates corresponding to (for 
the review of different models of space-time evolution of 
wave packages, see Ref. HH) 



l coh = (0.3 - 0.4 fm) • Ph /(GeV/c). 



(8) 



In the case of pions, it is checked to a large extent by 
the measurement of the A and Q 2 dependence of trans- 
parency in the A(e,e'n) reaction [36j . The condition of 
Eq. (J7J) is obviously much less stringent than for the case 
of CT phenomena with nuclei. It corresponds to require- 
ment that final momenta of the two hadrons with large 
Pt are 



p c >3^4GeV/c, p d > 3 4 GeV/c, 



(9) 



which is satisfied for the incident nucleons with pjq > 8 
GeV/c. 

Spins and parities of N and A are \ + and | + , so that 
the state h should be in 1 + . If it is a p-wave pion- or rho- 
like state, it is allowed. In general, h could be expressed 
by the Fock-space expansion, \h >= \qq > +\qqqq > H — • 
However, higher Fock states are suppressed because more 
hard gluon exchanges are needed to form the final pion 
and nucleon with large momenta. Therefore, the inter- 
mediate state should be considered as a qq state. Gluonic 
states such as gg are excluded as the intermediate state 
h for N — > A because the isospin state of h should be 
one by considering isospins of N and A are \ and |, 
respectively. They are present for the B = N, though 
the data indicate [32| that contribution of diagrams with 
gluon exchanges is much smaller than that with quark 
exchanges. 

If the final baryon (B) is a nucleon, there are ad- 
ditional contributions in principle. The same discus- 
sions on isospin, spin, and parity suggest that the in- 
termediate state should be in + and 1 + states. The 
1+ state is the same as the N — > A case. There is 



an additional intermediate state with the vacuum quan- 
tum number + . Therefore, the Fock-space expansion is 
\h >= \gg > +\qq > +\qqqq >+■•■. We do not take 
into account the + intermediate state in our current 
work with the following reasons. Lattice QCD calcula- 
tions indicate that the mass of the lightest + glucball 
is estimated about 1700 MeV (3^- Even if they were in- 
cluded in the cross-section estimate, their contributions 
are not significant in the kinematical region of \t\ <^ m 2 N 
because of their large masses. The /o(600) (or er) and 
/o(980) mesons are not S- wave qq states, so that their 
contributions are suppressed in the limit t' — > oo with 
t'/s'=const as follows. We have a mcson-nuclcon scat- 
tering amplitude which is proportional to the integral of 
the initial and final meson light-cone waves function over 
quark momenta: J 4>m (z,p±)dzd 2 p± , where z and p± arc 
the lightcone momentum fraction and transverse momen- 
tum of a quark, respectively. However, wave functions 
of mesons which are in the P-wave in the nonrelativis- 
tic limit have the property 4>m(\ — z,p±) = —4>m{z,Pa_) 
[40| . so that the integral vanishes since the final-state 
meson wave function is symmetric. If the fo mesons are 
tetra-quark (qqqq) or meson- molecule states 4l|, their 



contributions are suppressed due to more hard gluon ex- 
changes. In any case, it is possible to take channels, such 
as p + p — > p + 7r + + n, where the vacuum-like exchange 
does not contribute. In this way, we consider only the qq 
state as the intermediate state: 



h = qq. 



(10) 



We would like to emphasize here that the practical 
problem we face is that though it is possible to write 
the amplitude as a sum of the terms representing con- 
volutions of nucleon (transitional) GPDs with the hard 
blocks H of the qq-3q scattering, these blocks are not 
known so far and hence generic factorized expression is 
not useful for numerical calculations. Accordingly, we use 
an observation from the current studies of GPDs that a 
distribution in the ERBL region over z is to a very good 
accuracy symmetric and very close to that in the wave 
functions of 7r and p mesons which are close to the asymp- 
totic form starting with virtualities of a few GeV 2 [l(|. 
This will allow us to use the information about elemen- 
tary two body reactions. 

We explain our description of the nucleon and iV — > A 
transition amplitudes ^£^^hB in terms of the nuclconic 
and transition GPDs in Sec. IIV1 and the amplitude 
^hN->izN is explained in Sec. [Vj 



III. GENERALIZED PARTON DISTRIBUTIONS 
FOR NUCLEON AND N -> A TRANSITION 



Generalized parton distributions for nucleon 
and TV — > A transition 



Generalized parton distributions are introduced in the 
N — > B amplitude for the nucleon and N — > A for 
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transitions with neutral exchange in i-channel (p — > p, 
p — ► A + ) with other transitions related to these ones by 
the isospin symmetry. The N + N — ► N + ir + B cross sec- 
tion is factorized into two amplitude in Eq. (|5|). In this 
section, the first term M^^hB is explained. The inter- 
mediate state h should be considered as a qq state for a 
reaction with large- momentum transfer \t'\, \u'\ 3> m N . 
It is important to note that extraction of such a qq state 
with a light-cone separation is related to the GPDs in the 
ERBL kinematical region [HI, Therefore, the ampli- 
tude ^#jv^/jb is simply given by AAn-*b times the factor 
describing a light-like separated qq pair. 

The GPDs for the nucleon is given by off-forward ma- 
trix elements of quark and gluon operators with a light- 
cone separation between nucleonic states [HI EJ. The 
matrix element of vector current is associated with un- 
polarizcd GPDs for quarks in the nucleon: 

M V N = J ^e lXx (N, Pe |^(-An/2)#(An/2)| N,p a ) 

= I N i> N ( Pe ) [H(x, £,t)t + E(x, £, t) ^ a ip N (p a ), 



2m n 




t = A' 



FIG. 6: Kinematics for the N — > B transition and generalized 
parton distributions. 



The generalized scaling variable x and a skewdness pa- 
rameter £ arc defined in terms of light-cone momentum 
k+ (= (k° + k 3 )/y / 2) and P+ as 



(11) 



xP + , A~ 



f 



-n ■ A, 



(16) 



where 7jv is the isospin factor for the nucleon given by 
I N =< l/2||T||l/2 > (\M N : lm\\M' N )/V2 H|: 



In = 1, v2 for p 



p, n, 



(12) 



where k is the average momentum of the quarks, A M is 
the four momentum transfer from the initial nucleon (a) 
to the final nucleon (e) (or A in the N — » A transition), 
and t is its squared quantity: 



a aP is defined by a a/3 = (i/ 2 )[7 Q , and is a light- 
cone null vector which satisfies 



Pe 



Pa, t=(j> e -p a f 



(17) 



n 2 = 0, n-P=l for P= -(p„+p e ). (13) 

The vector is chosen so as to have n + = n_i_ = 0, 
where is the transverse vector and is defined by 
= (n° ± n 3 )/^/2. Then, is simply denoted as n M 
in the following discussions 0, [44| : 



V2P+' °' °' 2P^ 



(14) 



The momentum P is the average over the initial and 
final nucleon momenta. The integral variable A is defined 
by A = P ■ z ~ P+z~. In Eq. (JTTJ) , ^(An/2) is the 
quark field. The functions H (x, £, t) and S(x, £, i) are the 
unpolarizcd GPDs for the nucleon. The Dirac spinor for 
the nucleon is denoted as tpNip)- Quark-spin dependent 
GPDs are also defined in the similar way for the nucleon: 



M 



N 



dX 
2^ 



e lXx (N,p e \^(-\n/2)M(Xn/2)\N, Pa ) 



lNtp N (Pe)[H(x,£,t) ^75 +E(x,£,t) U n Al5 ] 1pN(Pa), 



2th n 



(15) 



where H (x, £, t) and E(x, £, t) are the polarized GPDs for 
the nucleoli. 



A schematic diagram is illustrated in Fig. [5] for under- 
standing the kinematics of the GPDs. The variable x 
indicates the momentum fraction carried by a quark in 
the nucleon as usually used in parton distribution func- 
tions (PDFs). The skewdness parameter £ or the mo- 
mentum A indicates the momentum transfer from the 
initial nucleon to the final one or the one between the 
quarks. In the forward limit (A — ► 0), the nucleonic 
GPDs H(x,£,t) and H(x,£,t) become usual PDFs for 
the nucleon as shown in Sec. IIII Bl 

Next, we explain how to express the transition ampli- 
tude A^tv^a in terms of the transition GPDs for N — > A. 
First, the N — ► A spin-independent quark distributions 
arc defined by a product of quark fields at a light-cone 
separation between nucleon and A states [Toj : 



M 



JV->A 



dX 
2^ 



e tXx (A, Pe |V(-An/2)^(An/2)| N,p a ) 



lANi>l(pe) [H M (x, f , t)K"n v + H E (x, £, t)¥Ll v n v 



H c (x, £, t)K% v n v ] ip N (p a 



(18) 



where Iaat is the N — > A transition isospin factor. The 
terms K.zf v , IC^ V , and K,^ v indicate covariants associated 
with magnetic dipole, electric quadrupole, and Coulomb 



G 




(a) - 1 < x < 



(b) 



(c) t, < x < 1 



FIG. 7: Three regions of x for the GPDs: (a) —1 < x < — £, emission and reabsorption of an antiquark, (b) — £ < x < £, 
emissions of a quark and an antiquark, (c) £ < x < 1, emission and reabsorption of a quark. 



quadrupolc transitions (451 ] : 



K 



M 



3(m A + m N ) 



ICf, 



-K 



2m N [(m A + tojv) 2 - t] ' 
M 6(m A + m N ) 



c 



^ m N Z(t) 
3(m A + mjv) 



A^P„ - A • PA„) 7 5 , (19) 



where m A is the A mass, and Z(t) is defined by 



Z(t) = [(m A + m N ) 2 - t][(m A - m N f - t}. (20) 



Our conventions for the antisymmetric tensor Sp^xa- and 
5 are taken as e i23 = +1 and j 5 = 170717273 OH 
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EH, The Rarita-Schwinger spinor [47J for the A is 
denoted as tp A (p)- 

The isospin factor Jatv in Eq. (JTSJ) comes from the 
transition isospin. The N — > A transition isospin is de- 
fine by [H 



|m 4 



T 



\m n 



^ 2\2 



( -Mat : lm 



(21) 



where e m is the spherical unit vector, < 3/2||T||l/2 >= 
2 is a reduced matrix element, and < 1/2 Mn : 
1to|3/2Ma > is a Clebsch-Gordan coefficient. The 
isospin factor becomes 



/a* = 1, yf, forp^A++,A+, A , A", (22) 



respectively. 

In Eq. (fl~8|) . the functions ifjvr(x,£, iJg;(x, £, f), and 
Hc(%,£,t) are the unpolarized GPDs for the transition 
JV — > A. Quark-spin dependent GPDs are defined in the 



same way [T(| : 

M N ^ A = j {A lPe \lp(-Xn/2)^(Xn/2)\N,p a ) 



Ian ^ A (p e ) 



A M (n ■ A) 



Hi(z,£,f)7v + if 2 (x,£,i)-^ 



< jV 



-A, 



TOat 

P • An,,, - 2A, 



'A' 



^Ar(Pa), (23) 



where the functions ifj(x,£,t) (i=l, 2, 3, or 4) [48| are 
the polarized GPDs for the N — > A transition. 

Partonic interpretation of the GPDs is different de- 
pending on x regions as illustrated in Fig. [7J The x 
variable is divided into three regions, — 1 < x < — £, 
— £ < x < £, and £ < a; < 1, and the GPDs are inter- 
preted as follows in each x region [ll[ : 

(a) "Antiquark distribution" 

-1< x < -£ (x + £ < 0,x - £ < 0) 

• Emission of antiquark with momentum fraction 
£-x 

• Absorption of antiquark with momentum fraction 

Both momentum fractions x + £ and x — £ are neg- 
ative, so that an emission (absorption) of a quark 
with the negative momentum fraction x + £ (x — £) 
corresponds to an absorption (emission) of an anti- 
quark with the positive momentum fraction — x — £ 
(-x + £). 

(b) "Meson (quark-antiquark) distribution amplitude" 



-£ < x < £ (x + £ > 0, x - £ < 0) 

• Emission of quark with momentum fraction x - 



-£ 



• Emission of antiquark with momentum fraction 

£-* 

The momentum fraction x — £ is negative, so that 
an absorption of a quark with the negative momen- 
tum fraction x — £ corresponds to an emission of 
an antiquark with the positive momentum fraction 
-x + £. 
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(c) "Quark distribution" 

£ < x < 1 (.t + £ > 0, x - £ > 0) 

• Emission of quark with momentum fraction x + £ 

• Absorption of quark with momentum fraction 



the quark contribution to the nucleon spin, possibly as 
well as the gluon one, is small, so that orbital angular mo- 
menta play an important role as the origin of the nucleon 
spin. These facts indicate that the GPDs are important 
for understanding basic properties of the nucleon from 
low to high energies. 



The regions (a) and (c) are called as DGLAP 
(Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) regions, 
and (b) is called the ERBL region [I?) depending on 
evolution types. 

If the reaction N — > h(qq) + A in Fig. [4] is dominated 
by the kincmatical region — £ < x < £, the cross sec- 
tion ofp + p^p + ir + A should be expressed by the 
transition GPDs for N — ► A. They correspond to the 
distributions in the kincmatical region (b) in Fig. [7J The 
GPDs for the nucleon are calculated, for example, in a 
chiral soliton model, and results indicated that the GPDs 
are dominated by a meson-pole contribution at small mo- 
mentum transfer and that the distributions are located 
dominantly in the region — £ < x < £ [49]. This supports 
our treatment of the blok N — ► h(qq) + A, in the calcu- 
lation of the process iV + iV^A^ + Tr + A, in terms of 
the transition GPDs. 



B. Properties of generalized parton distributions 
for nucleon 



Here we introduce basic properties of the GPDs for the 
nucleon. In the forward limit (£ = t = 0), the GPDs H 
and H for quarks are given by 

H(x, 0, 0) = q(x), H{x, 0, 0) = Aq{x). (24) 

where q(x) and Aq{x) are unpolarized and polarized par- 
ton distribution functions in the nucleon. Next, their first 
moments arc the form factors of the nucleon: 



dxH(x,£,t) = F^t), f dxE(x,£,t) = F 2 (t), 
-l J-i 

dxH(x,Z,t)=g A (t), f dxE{x,Z,t)=g P {t), (25) 
-l J-i 



where Fi(t), F2(t), gA(t), and gp(t) are Dirac, Pauli, 
axial, and pseudoscalar form factors of the nucleon. A 
second moment gives a quark orbital-angular-momentum 
contribution (L q ) to the nucleon spin: 

Jq = 2 / dXX ^ ^ ^' i = °) + E ( x > 1 = °)] 



= 2 A 9 + L 9> 



(26) 



where Aq is the quark-spin contribution and J q is the to- 
tal angular-momentum of quarks. It indicates that nuclc- 
onic GPD studies are important for clarifying the origin 
of the nucleon spin. As mentioned in the introduction, 



Properties of generalized parton distributions 
for N — > A transition 



In this section, basic properties of the N — > A transi- 
tion GPDs are explained. First, the first moment of the 
GPDs corresponds to the N — > A transition form factors 

MM-- 



J dxHi{x,£,t) =2Gi(t), i = M,E,C 



(27) 



where C?m GE(t), and Gc(t) are the magnetic dipolc, 
electric quadrupole, and Coulomb quadrupole transition 
form factors, respectively The N — > A transition is dom- 
inated by the Ml transition; however, there are also small 
E2 and C2 contributions 0, [24[ ■ Experimental investi- 
gations have been done in electron scattering, particu- 
larly in the processes such as (e,e'n) [22l [24 |50||. The 
E2 and C2 transitions are especially important for inves- 
tigating a dynamical aspect of baryon structure because 
they are related to the deformation of the nucleon and 
A. It is interesting to note that information on such a 
deformation should be also contained in the GPDs. The 
studies of the TV — ► A transition GPDs should be valuable 
for such a new aspect of low-energy hadron structure. 

The first moments of Hi(x, £, t) (i=l, 2, 3, 4) are axial 
form factors for the N — > A transition: 



1 r 1 

dxH 1 (x,a i t)=2C£(t), / dxH 2 (x,t,t) =2C£{t), 

-l J-i 

1 r 1 

dxH 3 {x,Z,t)=2C£(t), / dxH i {x,i,t)=2Ci{t), 

-l J-i 

(28) 



where Cf(t)(i=3, 4, 5, or 6) are the axial form fac- 
tors in Ref. j51j. These axial form factors were studied 
in parity- violating electron scattering (52l |. These Adler 
form factors are estimated in various theoretical models 
[52l ]. For example, a recent estimation by the chiral per- 
turbation theory indicates that C$ and Cf vanish in the 
leading order and they become finite due to higher-order 
corrections: \C£(t)\, \C£(t)\ < 1 [E|. 

Next, the TV — > A GPD studies are also valuable 
for understanding the origin of the nucleon spin, espe- 
cially orbital-angular-momentum effects. As an exam- 
ple, possible connections to the nucleon spin are shown 
by taking the large N c limit. In the chiral-soliton model, 
the N — > A transition GPDs are related to the nucle- 
onic GPDs because both baryons are different rotational 
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states of the same soliton fic| : 

H M (x, £, t) ~ [E u (x, £, t) - E d (x, i, *)] , 
Hx{x, C, t) =VS [S u (x, £ t) - H d (x, £, t) 

H 2 (x, £, t) =^ [£%t, t) - £ d (x, £, t) 
H E {x,Z,t) =H c {x,Z > t)=H 3 {x ) Z,t) 



(29) 



in the leading order of 1/N C expansion. It is especially 
noteworthy that the second moment of H M (x, £, t) is re- 
lated to the isovector part of the angular momentum of 
the nucleon carried by quarks as 



lim 

t— t0,JV o — >oo 



dx xHm(x, £, t) 

-l 

-j=[2(r-J d )-M? + Mi], (30) 



where 



£ [</(a:) + g(x)] , (31) 



be included in calculating the cross section. The calcula- 
tions are done in the same way as mesonic contributions 
to antiquark distributions 

Using the spin summation of the Dirac spinors, 
J2x N ^n{p)^n{p) =i>~ m N and Eq. (JTTJ for M N , and 
then calculating traces, we express the vector part of Eq. 
(f3"2"|) in terms of nucleonic GPDs: 



E \mw = i$ 



8(l-e){H(x,U)} 2 



+ WeH(x, t)E(x, Z,t)-—(1 + 2 {E(x, e, t)Y 



(33) 



In the same way, the axial term is calculated as 



E \m a n \ 2 = i. 



8(i-^ 2 ){i?(x,e,t)} 5 



18Z 2 H(x,Z,t)E(x,S,t) 



2te 



{E{x^,t)Y 



N 



(34) 



These expressions are used for calculating the cross sec- 
tion for N + N^N + n + N. 



is a second moment of quark and antiquark distribution. 
The J" and J d indicate total angular momenta carried by 
u and d quarks, respectively. They are equal to the sum 
of the spin and orbital angular momentum carried by the 
quark, as shown in Eq. (|2l)|) . Equation (|3"0|) suggests that 
the N — > A GPDs should be important for investigating 
the isovector part the orbital angular momentum, L u — 
L d , in the studies of nucleon spin. 



IV. N -> B TRANSITIONS IN TERMS OF 
GENERALIZED PARTON DISTRIBUTIONS 

A. Generalized parton distributions for the nucleon 

For calculating the unpolarized cross section of N + 
N — > N + 7r + N, spin summations are taken for the nu- 
cleons. In this section, the factor describing transition 
from the initial nucleon (a) to the final one (e) is calcu- 
lated. We discussed the amplitudes and M^, which 
are described by spin independent and dependent gen- 
eralized quark distributions. The absolute- value squared 
of the transition amplitude is given by 



E i^i 2 = E (\ M 



V |2 



\M N f) 



(32) 



It should be noted that there is no interference term if 
the spin summations are taken over Xn (= X a ) and Aap 
(= A e ). The summation of the vector and axial contri- 
butions indicate that all the associated qq states should 



B. Generalized parton distributions for the N — > A 
transition 



The N — ► A transition part is calculated in the same 
way by taking spin summations for the nucleon and the 
A. The N — > A transition was described by the matrix 
elements M-Yf^A and M-n->a with the transition GPDs 
in Sec. IIII Al Then, absolute-value squared of the tran- 
sition amplitude is 



E 

Aw , Xa 



\M 



N^A 



E 



\M 



A jy , A/ 



(35) 

Because the spin summations are taken over A at and Aa, 
there is no interference term. 

In our numerical estimation for the cross section, we 
use the following approximations. First, the electric 
quadrupole contributions He and He should be small 
as suggested in theoretical and experimental studies in 
the N — > A transition form factors [13, HJ- Therefore, 
it is appropriate to assume He(x, £, t)=Hc(x, £, t)=0 for 
our first estimation of the cross section if such small ef- 
fects can be neglected. This is also satisfied in the chiral- 
soliton model as shown in Eq. (|2"9"|) . Then, the vector 
amplitude is expressed only by the magnetic term: 

■M-N^A = lNA'<pA(Pe)HM(x,Z,t)lCf^n V 1pN(Pa)- (36) 

Taking the square of its absolute value and performing 
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spin summations, wc obtain 

E \M^ A \ 2 = -iA{H M (x,c,t)y 



Ajv ,A, 

x Tr 



^^iPe)^l{Pe)^ u n"^i) N (p a )i) N {p a )lC^ 

. Aa Ajv 



(37) 



Spin summations of the Rarita-Schwingcr spinors are 
given by [13] 



forms are 



Ci (£, t) = ^-5- (1 - 2 [(m A + - t] , (43) 



Ca(£,t) 



£ 2 [(m A + mw) 2 - <] 



x [ - 4m 2 A t + (m A - m|r + i) 2 ] , (44) 
8£ 



[(mA + mjv) 2 - t] 



x [4m|£ + (1 - - mf, + t)] . (45) 

In this way, the iV — > A transition part of the cross 
section can be calculated using the transition GPDs. 



3?7lA 



3m A ' ' 



(38) 



hN -> ttN SCATTERING AT LARGE 
MOMENTUM TRANSFER 



where m A is the A mass. We denote the trace part of 
Eq. dnZD as C M (£,*): 

E l^-Al 2 =^ A CM(e,t){^M(a:,^t)} 2 , (39) 

Atv,Aa 

and it is calculated by using Eq. (j3"5)) as 



771AT + 771A 



Cm(£,*) - 3 

_mAf{(mAr + ?tia) — t\ 
x [t- (tun + m A ) 2 ] 

x [t(l - C 2 ) + 2£(m| - m 2 w ) + 2^ 2 (mi + mfv)] . (40) 



Similarly, the axial-vector part is given by (neglecting 
the GPDs H 3 (x,£,t) and H^x^.t) as suggested, for ex- 
ample, by the chiral-soliton model [l(J El and also by 
the chiral perturbation theory of the form factors [HI): 



A /J (n • A) 



N 



i>N(Pa). (41) 



Expressing these amplitudes in terms of the GPDs in 
Eqs. ((2)]) and we obtain 



E ^ 



A \2 _ j 2 
AT— > A I — J JV A 



Ci(e,t){ffi(a;,^t)} 2 



+ C 2 (e,t){ J ff 2 (x,^,t)} 2 



(42) 



where Ci(£,t) are coefficient terms and their explicit 



A. Intermediate qq and mesons 



As explained in Sec lIII M the GPDs are separated into 
three x regions. If the momentum transfer \t\ is small, 
the GPDs distributed dominantly in the — £ < x < £, 
for example, as suggested in Ref. [4!|. The GPDs in 
the kinematical region, — £ < x < £, are the distribution 
amplitudes with the quark and antiquark emissions with 
momentum fractions x + £ and —a; + £, respectively. 

In general, all associated qq states should contribute to 
the cross section as we considered the vector- and axial- 
vector amplitudes. At small \t\ <C mfy, light-meson-likc 
qq pairs contribute mostly to the cross section. The dom- 
inant ones are the qq state with the pion (p) quantum 
numbers in the axial-vector (vector) amplitude. Then, 
the qq pair could be considered as a pion or p state in 
a small-size configuration and the amplitude is given by 
the meson (M = tt or p) wave function <Pm{z) 1 and its 
normalization is given by the Brodsky-Lepagc relation 
[55l [50!] . However, our cross sections are independent of 
the functional form of this wave function in the region 
where interpolation of the GPDs by the lowest states is 
a good approximation (see Sec. IVI[) . The kinematical 
variable z is the light-cone momentum fraction and it is 
defined in the following way. The variables x\ and x 2 
indicate momentum fractions carried by a quark and an 
antiquark, respectively, in the nucleon or the N — > A 
transition. As indicated in Fig. [6l they are expressed by 
x and £ as 



xi = x + £, x 2 = -x + £. 



(46) 



The variable z is the momentum fraction carried by the 
quark to the total momentum of the quark and antiquark 
system: 



X\ + x 2 



(47) 
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In the previous section, the N — > B part was calculated 
in terms of the nucleonic and N — ► A transition GPDs, 
which are associated mainly with the qq emission. In the 
small momentum transfer region (|i| <C Tn%), the qq pair 
could be considered as a pion or p state in a small-size 
configuration, which next interacts with another proton 
(b) as shown in Fig. |U The probability that the qq 
pair forms a meson M is described by the meson wave 
function in the light-cone coordinates. However, it should 
be noted that the GPDs at the meson pole effectively 
contain the wave function of the meson in the point like 
configuration {4>m{z)) as indicated, for example, in Ref. 
[l(| as the pion-pole contributions. This factor enters 
also in the meson- nucleon scattering amplitude, so that 
this factor should be removed in matching of the GPD 
description with description in terms of scattering off a 
meson: 



n iK,p', 0) 

A 



,90° 



M(E M ,0,p) 



p(E p ,0,-p) 



p(E' p ,-p',0) 
FIG. 8: Mp — > -np elastic scattering at 6 c . m . — 90° 



M 



Mp — >7ip j 



(48) 



where {qq)M is a qq state with the quantum numbers of 
the meson M. The critical point here is the observation 
of several GPD analyses that z distribution of qq pairs 
in the GPDs is close to that in the pion and rho mesons. 
The off-shell nature of the qq-N scattering amplitude can 
be best seen from s' = —t'—u'+2m N +m 2 .+t as compared 



2m 2 N 



ml + m 2 M , 



to the on-shell case of s' = — t' — u 
where tum is the meson mass. However, in the limit Eq. 
d3j) , the cross section is a function of the ratio of large 
variables and so correction goes to zero. 

For the spin-dependent GPDs in the axial amplitudes 
Mff and M^^ A7 the factors are approximated by the 
pion wave function as long as the qq pair has the same 
spin structure of the pion. For spin-independent GPDs 
in the vector amplitudes and M^_ tA . l the factors 
are related to the p meson. In this way, the TV + N — > 
N + ir + B reaction can be calculated by the N — > B 
transition part and the MN — > irN scattering process: 

■MnN-^NttB = ^N->h(qq)B ' -^h(qq)N-*TTN , 

M MN-urN /4>m{z), (49) 

at small momentum transfer <C mjy). 

It is worth noting here that the discussed process is 
directly sensitive to the dependence of the GPDs on £. 
Dependence on x is integrated over. This is similar to the 
situation in most of the other exclusive hard processes, 
for example the reaction 7* + p — > ir + + n. 

Now, the remaining part is the description of the M + 
TV — > 7r + N scattering. 



B. Parametrization of cross section for 
meson-nucleon elastic scattering 

In order to calculate the cross section of meson-nuclcon 
elastic scattering Mp — > irp, where p indicates the pro- 
ton, we try to obtain useful parametrizations of the cross 
sections by using experimental measurements. The cross 



section is first expressed in terms of parameters which 
are then determined by fitting np — > np and np — ► pp 
scattering data. The momentum-transfer dependence of 
the cross section is typically parametrized as [32j 



da 



Mp — >7rp 



(S>,t>) 



dt> 



oc a + c(t' — t' o y 



(50) 



where a and c are parameters, and t' is defined by the t' 
value at the scattering angle 9 C m = 90° in the center-of- 
mass (cm.) system. The Mandelstam variables s' and t' 
are given in Eq. (J5J). Moreover, this cross section scales 
with the cm. energy squared (s') as 



d(J Mp^np(s' ,t') 1 



dt' 



Mp 



(t'/s') 



for 



00 at fixed — , 
s' 



(51) 
Here, 



according to the counting rule [3l|, |38I |. 

fMp^Trp(t'/s') is a function of t'/s'. The factor n is the 
total number of all interacting elementary fields, and it is 
basically the sum of valence-quark numbers in the initial 
and final states. In the meson-nucleon elastic scattering, 
n = 10 and so the energy dependence is 1/s' 8 . 

The factor t' Q is defined in the cm. frame, in which 
our momentum assignment (E,pt,pl) is shown in Fig. 
[8] Here, pr and pl indicate transverse and longitudi- 
nal momenta, respectively. From the energy-momentum 
conservation, the energies and momenta are given for the 
mesons and protons as 



E p = 



{ml 



+ {ml 



2vs' 
2 



L M) 



2V? 



Ei 



(ml - ml) 



2vs' 
s' + {m 2 p - ml) 



P = 



^Js' 2 - 2s' {ml + m 2 M ) + (mjj 



4) 2 



s' 2 — 2s'{m 2 + m 2 ) + {m 2 — m 2 )' 



2\/7 



(52) 
(53) 

(54) 

(55) 
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TABLE I: Determined parameters in Eq. ([57)l . 





Process 




a (barn ■ (GeV) 2n " 6 ) 


Parameters 
c (barn ■ (GeV) 2n " 10 ) 


n - 2 (Ref. [32]) 


x7d-o.f. 


7T 




P 


2.08 ±0.10 


1.11 ±0.30 


6.7 ±0.2 


0.96 


7T 


~p — » TT~ 


"P 


8.94 ± 0.42 


10.46 ± 1.55 


7.5 ±0.3 


0.36 


7T 


+ P -> 


P 


189 ± 14 


60.0 


8.3 ±0.5 


0.21 


7T 


~P P" 


P 


218 ± 12 


200.0 


8.7 ± 1.0 


0.80 



3.5 



■ 71 + p— > 71 +p 

• Jt~+p— » 7t~ + p 




O 



□ 7l + +p— > p + +p 












5 } V5- U --' 





■t' (GeV 2 ) 



4.5 5 

-t' (GeV 2 ) 



FIG. 9: Data of n p — > n p and 7r~p — + n~p cross sections 
from the BNL-E838 experiment [12] and our parametrization 
curves of Eq. (|57|l with the optimized parameters. The solid 
and dashed curves are for the ir + p — > ty + p and 7r~p — > 7r~p 
cross sections, respectively. 



FIG. 10: Data of ix + p — > p + p and 7r~p — > p~p cross sections 
from the BNL-E838 experiment [32| and our parametrization 
curves of Eq. ([570 with the optimized parameters. The solid 
and dashed curves are for the ir + p — > p + p and 7r~p — > p~p 
cross sections, respectively. 



Therefore, t;q becomes 

s' 2 - s'(2m 2 + m 2 M + ml) + (to 2 - m 2 M )(m 2 p - ml) 



2s' 



(56) 



In this way, the differential cross section is parametrized 
in terms of the variables s' and t' as 



da 



Mp — >7rp 



(s',t>) 



1 



(if 



„/n-2 



[a + c{t' -t' a f 



(57) 



Experimental data are taken from the Alternating Gra- 
dient Synchrotron (AGS) experiment E838 [13] with the 
incident pion-beam 5.9 GeV at the Brookhavcn National 
Laboratory (BNL). They observed large-angle scattering 
cross sections for 7r + p — > 7r + p, n~p — > n~p, n + p — > p + p, 
and 7r~p — > p~V - By analyzing the scaling behavior be- 
tween E755 [53] and E838 data, the parameter n — 2 was 
determined as 6.7, 7.5, 8.3, and 8.7 in Ref. as listed 
in Table [J 

We use Eq. (|5"T|) for describing the cross sections by the 
exchange of initial and final states (Mp «-> np). The pa- 
rameters a and c are determined by fitting the E838 data. 
As shown in Figs. [91 and [TOl there are 16, 21, 6, 8 data 
points for n + p — ► 7r + p, n~p — ► 7r~p, 7r + p — > p + p, and 
7r~p — > p~p, respectively. In Table [H determined param- 
eters are shown and x 2 values per degrees of freedom are 



X 2 /d.o.f.=0.96, 0.36, 0.21, and 0.80. The two-parameter 
fit for the it p — > p p was difficult due to a small num- 
ber of data. Therefore, the parameter c is fixed at 60 
and 200 for 7r + p — > p + p and 7r~p — > p~p, respectively, so 
as to have similar i-dependency to the 7r ± p — > 7r ± p cross 
sections in Fig. [5] The fits are successful in reproducing 
the data except for some deviations at t ~ — 5.5 GeV 2 in 

Fig. [Hand at t 4.3 and -5.0 GeV 2 for tt^p -> p~p 

in 7r~p — > p~p of Fig. [TO] 



VI. NUMERICAL RESULTS FOR CROSS 
SECTIONS 

A. Expression for N±N^N±7r±B cross section 

In this section we calculate the cross sections for the 
N + N — > N + it + B using the parametrizations for 
the M ± N — > 7r ± N cross sections, and a model for 
the nucleonic and N — > A transition GPDs. Here, the 
expression for the cross section is summarized. 

The MN — > irN cross section is given by 



da 



1 



MJV-»ttJV 



d 3 p c d 3 p d 



4 v / (Pb-PM) 2 -m 2 r p 2 f 2£ c (2tt)3 2i? d (2^)3 
^ I] \M MN ^ N \ 2 (27:) 4 5 4 (p b + p M - Pc - p d ), (58) 
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where ^2\ b indicate the average over the initial nucleon 
(b) spin. The general expression for the cross section is 
given in Eq. (JJ) for the process N + N — ► N + ir + B. The 
factor S is S — 1 in our numerical results for B = A , 
A ++ , or neutron, which is not identical to the proton. By 
assuming the factorization and mesonic qq dominance, 
the matrix-element part is separated into two processes 
as shown in Eq. ([15]) . Substituting Eq. (J49J) into Eq. 
(UJ) , we write the N + N^N + tt + B cross section by 
the N —> B amplitude, namely by the TV -> B GPDs, 
and the MN — ► ttN scattering cross section: 



da 



NN^NttB 



-EE 

A„ A„ 



V (Pb ■ Pm) 2 - m 2 N P 2 M d3 Pe 
V(Pa-p b ) 2 ~rn% 2E e (2nr 

i 

;\M N ^ B \ 2 da M N->-xN{s' ,t'). (59) 



We take the z axis as the momentum direction of the 
initial nucleon b, so that our momentum assignment for 
p a ,Pb, and p e is 

Pa = {m N , 0,0,0), 
p b = (E N ,0,0,p N ), 

p e = (E B ,p B sm9 e cos(/) e ,p B siii9 e siii<j) e ,p B cos9 e ), (60) 

where 9 e and <p e are the polar and azimuthal angles for 
the baryon B. The integral over the B momentum p e is 
written as 

d 3 p e = pldp e dtt e = p 2 e dp e d(cos 9 e )d4> e , (61) 

The differential cross section for NN — > NirB is ex- 
pressed with respect to the kinematical variables t, t' , y, 
4>d, and </> e , where the variable y is defined by 



V 



t + m 2 N + 2(m N E N - E B E N + p B PN cos 9 e ) 



(62) 

The variable 4>d is the azimuthal angle for the final nu- 
cleon (d). Then, integrating over y, tpd, and </> e , we finally 
obtain the following differential cross section: 



da 



NN^NttB 

dtdt' 



dy 



16 (27r) 2 win Pn 



(ys-t - m N ) 2 - Am 2 N t da M N^Ni.s' = ys,t') 
(s — 2m 2 N ) 2 — Am\ 



N 



dt> 



■\M 



N^B\ , 



(63) 



where y ma x and y m in are upper and lower bounds of the 
integral. From Eq. (|62|) . y m ax is obtained by taking 
cos# e = 1- The lower bound is estimated by applying 
a hard kinematical condition for t' and u': \t'\, \u'\ > 
Qq, where the hard scale is taken in the region Qq=2~3 



GeV 2 . The relation s' + 1' + u' = t + 2m 



N 



i 2 leads to 



Vmin — 



Qo 



2m% - , -t'>Ql (64) 



where \t\, m 2 -C m 2 N is assumed. 

We specify all used kinematical variables and con- 
stants. The initial proton energy is given by 



E N = m N + T 1 



N- 



(65) 



where the kinetic energy of the proton beam (Tjv) is, for 
example, 30 GeV at the initial stage of J-PARC and is 
expected to be 50 GeV at the later stage. Then, the cm. 
energy squared is calculated as 



s = 2m 



N 



2mjqE 



AT) 



(66) 



The momentum transfer squared t is defined in Eq. ([2]) 
and we present the cross section as a function of this 
variable in Eq. (|6"5|) . The value of t determines the energy 
and momentum of the final baryon B: 



Er 



t 



2mjv 



Pb 



' B- 



(67) 



in the rest frame of the nucleon (a) . The polar angle of 
Pb is given by 



COSfA, = 



ys-t - m 2 N + 2{E B - m N )E N 
2pb Pn 



(68) 



In calculating the skewdness parameter £, we should 
be careful about the reaction kinematics. We consider a 
reaction by taking z axis along the incident proton beam, 
and the GPDs of a target hadron at rest are studied. In 
the cm. frame, the target momentum along the z axis is 
negative. Therefore, the appropriate choice of the vector 
n M should be n/f instead of the usual in many GPD 
articles: 



1 1 

-, 0, 0, 



V2P-' ' ' V2P- 
The skewdness parameter £ becomes 



1 m N - E B +Pb cos 9 e 

c = 71 ■ Za = 

2 7727V + Eb — Pb cos 6 e ' 



(69) 



(70) 



by using A = p e — p a and Eq. (j60|) . Substituting Eq. 
into this expression, we obtain 

ys-t- m 2 N + 2(E B - m N )(E N - p N ) 



ys - t - m 2 N + 2{E B - m N )E N - 2(E B + m N )p 



N 



(71) 



The kinematical minimum £ m i„ is calculated by substi- 
tuting y m i n into the above equation, whereas the maxi- 
mum ^max is obtained by Eq. (|70p with cos0 e = 1. 

Here, it should be noted that the energy and momen- 
tum of B depend on the variable t as explicitly shown in 
Eq. (|67[) . The relation between £ and y is obtained from 
Eqs. (|62 l) and (|70 | : 



1 

y = - 

s 



t + m 2 N -2{E N -p N )E B 
+ 2mjv£'A r 



2 mN PN 



(72) 
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which will be used in changing the variable y for £ when 
presenting the numerical results for the cross section in 
Sec. rVICl 

The elastic scattering factor daMN^nN(s',t')/dt' is 
calculated by using the result in Sec. IV B| in particular 
Eq. ([FT]) together with the parameter values in Table Q] 
The last ingredient is the matrix elements Mjv— >b which 
are expressed by the GPDs. To be specific, we use Eqs. 



\M 



N 



The 



([55 ]) . and ([53 ]) for calculating ^ Ajv . 
TV — > A matrix element J2\ \ \Mn^a\ 2 is calculated 
by using Eqs. (09]), {gjj| g3J, (gSJ), ©, and gSJ. 

As we explain in the following subsection, in our numeri- 
cal analysis we estimate the meson-pole contributions to 
the cross sections. The meson-pole contributions to the 
GPDs contain the meson wave function (f>M(z), which 
is divided out from the cross-section expression in Eq. 
(f6"3")) . Therefore, a specific functional form of 4>m{z) is 
not necessary for our current numerical estimates. 



B. Contributions from the meson poles 

In this article we want to provide first estimate of the 
order of magnitude for the N + N — > iV + 7r + _B cross sec- 
tions in order to establish feasibility of investigating the 
GPDs at hadron experimental facilities. For this purpose, 
we employ a very simple description by meson poles. Our 
studies should be intended merely to provide a starting 
point for future theoretical and experimental studies. 

In the ERBL region, major contributions to the GPDs 
come from meson poles, particularly the pion and p- 
mcson poles, at small \t\. For numerical evaluations of 
the N + N^B + tt + N cross sections, we show explicit 
expressions for these pole contributions to the nuclconic 
and N — > A transition GPDs. 

First, the pion is a pseudo-scalar particle, so that it is 
related to the axial part of the nucleonic matrix element. 
The pion pole contributes only to the nucleonic GPD 
E(x,£,t) as [MIH 



^9A m 2 N 



V3 



1 - t/ml ml fa £ 



fa(z)0(Z-\x\), (73) 



where gA is the axial charge of the nucleon and its value is 
gA = 1.2695 (58|. Here, E indicates the isovector combi- 
nation E u — E d . In the same way, p contributions can be 
estimated. Since the p is a vector particle, it contributes 
to the vector part of the matrix element. The p-pole 
terms are given by extending the expression in Ref . [59j : 



1 d y/3 



1-t/m* g P V2/ P £ 



4> p {z)e^-\ x \), 



^^^d^g^^^-N), (74) 

where the G\ and C 2 are two-types of p coupling con- 
stants with the nucleon, and they are given for the isovec- 



tor couplings, namely for the isovector p, by 

C\ e p — e n 1 C2 n p — n n 1.85295 
9 P 2 2 ' g p 2 • 2m N 2m N 



(75) 



where e p (e n ) is the charge of the proton (neutron) and 
k-p (k„) is the anomalous magnetic moment of the proton 
(neutron) (58j . There is no pion contribution to H (x, £, t) 
because the function H is chiral even, whereas the pion 
should contribute to the chiral odd. 

We may note that normalizations of the pion and p 
wave functions arc given by 



1 = 



1 /o 

dx9^-\x\)^—4^{z) 



dx6(£~- \x\) 



V2/ P £ 



(76) 



where the normalization factors are slightly different from 
the ones in Ref. [lol ]. 

As for the light-cone wave function of the meson M, 
we may take, for example, asymptotic wave function of 
pion and p mesons (55l . |56| : 



<j> % (z) = V3fnz(l - z), 
(f) p (z) = V6f p z(l - z). 



(77) 



The normalization of these wave function is uniquely 
fixed in QCD [H, [56[. For the pion wave function it 
is given by 



dZ(j)Tr(z) 



dz d 2 k± 
16tt 3 



2V3' 



(78) 



where / ff is the pion decay constant for 7r + — > p + fp. 
By noting the difference of \/2 from the value in Ref. 
[Ha, it is /tt=92.4 MeV. The p decay constant f p is 
determined by the decay p° — > e + e~, and it is given by 
ii 



3m p r p a_ e + e - 



8na z 



= 110.6 MeV. 



(79) 



where m p is the p mass, a is the fine structure constant, 
and r p o_> e + e - is the p° — > e + e~ decay width. 

Next, pion and p contributions are shown for the ./V — > 
A transition GPDs. The pion contributes to the axial 
part, particularly the GPD H2 S3, is given in the 
same way by 



%9A m 2 N 



(1- t/ml) ml 



Mz)0(Z-\x\). (80) 



The p-meson contribution to the vector part is expressed 
by the magnetic N — > A transition as 



G 



M 



V3 



I -t/ml V2f P t 
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where G* M is the transition magnetic moment [45|. Ex- 
perimental measurements indicate G\ [ rss 3 (50l ] , which is 
close to an expectation G* M — (/i p — /x n )/V2 = 3.32757 
[Io| , where and are magnetic moments for the pro- 
ton and the neutron, respectively. In our numerical anal- 
ysis, we use the value G* M = 3.3 for the transition mag- 
netic moment. 



C. Numerical results 

From measurements of the N + N^B + n + N cross 
sections, the GPDs should be extracted from the data. 
There is no such data at this stage, so that we need to 
provide an order of magnitude estimates of the cross sec- 
tions for experimental proposals at hadron experimental 
facilities, for example, J-PARC and GSI-FAIR [H, 
Since large contributions should come from pion- and p- 
pole terms, such pole contributions are evaluated in this 
work. However, one should note that the following nu- 
merical evaluations based on the meson-pole terms are 
just rough estimates of the cross sections. In particular, 
there exists an extra t dependence in the GPDs at large 
\t\ (like form factors at large \t\) in addition to the pole t- 
dependence, and it is neglected in this work. Conversely, 
the actual GPDs need to be extracted from future mea- 
surements once the data are taken. Then, they should 
be used for investigating the spin structure of the nu- 
cleon, especially orbital angular momentum effects, and 
electromagnetic properties of the N — ► A transition as 
explained in Sec. IIII Bl and Sec. IIII CI 

In the case when incident particles are unpolarized and 
polarization of the final particles is averaged over, we 
predict the same cross section for analogous reactions 
initiated by antiprotons. Interference of axial and vector 
terms survives if the target proton is transversely po- 
larized. Interestingly, the sign of the interference effect 
for incident protons and antiprotons depends on relative 
phase of the A^n-t^n and Ap^^^N amplitudes. In lead- 
ing QCD, diagrams for this process correspond to real 
amplitudes, so the relative phase could be either 1 or —1. 
So, in principle, by combining measurements with proton 
and antiproton beams, one can measure this phase. 

We show the numerical results for the following three 
cross sections 

- + A°, 

- + A++, 



(1) p + p — ► p + 71 

(2) p + p — > p + 71 

(3) p + p — > p + 7T+ + n, 

in Fig. QT] for the kinetic energy of the proton beam at 
30 and 50 GeV. We focus on these processes for several 
reasons: they are among the simplest to measure exper- 
imentally, they are expressed through the GPDs which 
are best studied theoretically, and they are expressed 
through elementary large angle meson-nucleon reactions 
which were studied experimentally. 

Also, in these processes the intermediate state with 
the vacuum quantum number + does not contribute as 



explained in Sec. HT1 which further simplifies the descrip- 
tion. The cross sections are calculated at the momentum 
transfer t = —0.3 GeV 2 , and the momentum cutoff in 
Eq. is taken as Qg=3 GeV 2 . The cross sections are 
decreasing function of — t', and they are of the order of 
\x barn/ GeV 2 . Hence it appears feasible to accumulate 
large statistics for these reactions over the range of t' for 
several incident energies and check the suggested mecha- 
nism of the reaction. Therefore it should be in principle 
possible to measure the GPDs in the N+N -> N+n++B 
reactions. The A-production cross sections are as large 
as the nucleonic one, which indicates that TV — > A transi- 
tion GPDs could be also studied as well as the nucleonic 
GPDs. At Tjv = 50 GeV, the magnitude of the cross sec- 
tions becomes slightly smaller; however, the larger beam 
energy is valuable for extending the kinematical region 
of t (and £ as shown later) and testing the mechanism of 
the reaction in particular by comparing the cross sections 
for the same s',t' and different Ejf. 

Next, the pion and p- meson contributions to the cross 
sections are explicitly shown for the processes p + p — > 
P + TT+ +n a,ndp + p ^ P + 1T+ + A in Fig. [12]at the 30 
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FIG. 11: Differential cross section as a function of t'. The 
incident proton-beam energy is 30 (50) GeV in the upper 
(lower) figure, and the momentum transfer is t = —0.3 GeV 2 . 
The solid, dotted, and dashed curves indicate the cross sec- 
tions for p + p — > p + 7r + + A , p + p^p + ir~+ A ++ , and 
p + p — > p + 7r + + n, respectively. 
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GeV beam energy and the momentum transfer t = —0.3 
GeV 2 . The results indicate that the contribution from 
the p-like vector GPD terms are smaller than the pion- 
like axial- vector terms in both process p+p — > p + 7r + + n 
and p + p^p + TT + + A°. Especially at large — t', the p 
contributions are much smaller than the pion ones. It is 
expected in general because cross sections should be sup- 
pressed for heavier-meson intermediate states. Therefore, 
it is important first to understand the axial- vector (quark 
spin independent) GPDs for estimating the magnitude of 
the cross sections. However, both the axial-vector and 
vector (spin dependent) GPDs should be calculated for 
estimating the cross sections. Once experimental data 
are obtained, they should be valuable for determining 
both spin independent and spin dependent GPDs for the 
nucleon and for the N — ► A transition. 

The cross-section results are presented in Figs. [TTIand 
[T2]for the hard momentum cutoff Q 2 =3 GeV 2 for t' and 
v! in Eq. (|64[) . We show the cutoff dependence of the 
cross sections in Fig. [T3] within the range 2 GeV 2 < Q 2 < 
4 GeV 2 . The cross sections depend strongly on the cutoff 
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FIG. 12: Differential cross section as a function of t' . The 
incident proton-beam energy is 30 GeV, and the momentum 
transfer is t = -0.3 GeV 2 . The upper (lower) figure indicates 
the cross section for the process p + p — > p + ir + + A° (p + p — > 
p + 7r + + n). The solid, dotted, and dashed curves indicate 
the cross sections for the total, axial-vector (n) contribution, 
vector (p) contribution, respectively. 
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FIG. 13: Cutoff (Qo) dependence of differential cross section. 
The incident proton-beam energy is 30 GeV, and momentum 
transfers are fixed at t = —0.3 GeV 2 and t' — —5 GeV 2 . 
The solid, dotted, dashed curves indicate the cross sections 
for p + p — * p + 7r + + A , p + p — > p + 7r~ + A ++ , and 
p + p — > p + 7r + + n, respectively. 
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FIG. 14: Differential cross sections are shown as a function 
of £ at the incident proton-beam energy 30 (or 50) GeV, t = 
—0.3 GeV 2 , and t' = —5 GeV 2 . The solid and dashed curves 
indicate the cross sections for p + p — > p + tt + + A and 
p+p — > p + tt + +n, respectively. The two upper (lower) curves 
indicate the cross sections at Tn = 30 (50) GeV. There is a 
restriction on the kinematical range, £ m i n < £ < £max, which 
is explained in the text. 



if it is taken low enough ~ 2 GeV 2 . As the momentum 
transfers \t'\ and \u'\ become smaller, soft contributions 
become larger. It will be important to study process as a 
function of —t' to determine minimal value of the cutoff 
which is necessary for using discussed hadronic processes 
for determination of the nucleonic and N — > A transition 
GPDs. 

In connection with the studies of GPDs, it is instruc- 
tive to present the cross sections as a function of £ with- 
out integrating them over the variable y. Using the re- 
lation between the variables y and £ in Eq. (fT2")l . we 
calculate the differential cross sections da / d^dtdt' . In 
Fig. [TJ] the cross sections are shown as a function 
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of £ at fixed * and t' (t = -0.3 GeV 2 and f = -5 
GeV 2 ). It should be noted that the kinematical range 
of £ is limited (£ m j„ < £ < £rna:c) as explained below Eq. 
(|71|) . The minimum £ m ;„ comes from the hard kinemat- 
ical condition for \t'\ and \u'\, and the minimum £ m i n is 
due to the minimum polar angle for the final baryon B 
(6 e = 0). Because of the larger masses for the A than 
the one for the nucleon, the £ range is more restricted: 
£ m ax = (m N - E B +PB)/(m N + E B ~Pb) for the A. At 
larger beam energy (Tjv = 50 GeV) , the smaller-^ region 
is probed by the reactions. Therefore, 50-GeV experi- 
ments arc valuable for extending the kinematical regions 
(in £ and also t) of the GPD measurements in addition to 
the 30-GeV experiments. Note in passing that the use of 
50 GeV proton beam would allow to produce secondary 
beams of mesons and perform precision measurements of 
meson-nucleon large angle scattering processes which en- 
ter in the analysis of the 2^3 processes and of interest 
on their own. 

Next, the same cross sections arc calculated at differ- 
ent momentum transfer points for t. At a given £, there 
is a restriction for the variable t at 6 e — 0, and it is 
given by {-t) min = 2£[m%/(l - £) - m 2 N /(l + £)]■ The 
cross sections are calculated at this point of t for each 
£ in Fig. [151 Then, the cross-section range is extended 
to larger £. Measuring these cross sections at future fa- 
cilities at J-PARC and GSI-FAIR, we should be able to 
obtain information on the GPDs in the ERBL region. 

In this article, we discussed processes N + N — > N + 
n + B (B=N or A) for investigating the GPDs at the 
hadron facilities. In addition to these GPD studies for 
the nucleon and A, strange-particle productions are also 
quite interesting, and they can be calculated along the 
same line. For example, the process p + p — » A+p + K + 
is likely to have a cross section comparable to that of 
the reactions we discussed in this work as the smaller 
KAN vertex is to some extent compensated by larger 
cross section of K + p — > K + p scattering. 
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FIG. 15: Differential cross sections are shown as a function 
of £ at the incident proton beam energy is 30 (or 50) GeV and 
t' = -5 GeV 2 . The notations are the same in Fig. [H] The 
momentum transfer t is taken at cos 8 e = 1 as a function of £. 



VII. SUMMARY 



We introduced a new class of hard 2^3 branching 
hadronic processes and investigated a possibility that the 
generalized parton distributions for the nucleon and the 
V — > A transition could be investigated in the reaction 
N + N — > N + it + B where B is the nucleon or A. 
First, we argued that the process is factorized into two 
blocks, N —> B and meson-nucleon (MN) scattering by 
imposing hard scattering condition for the meson-nucleon 
block. The V — > B part is expressed in terms of the nu- 
clconic and the N — » A transition GPDs, and the MN 
scattering part is parametrized so as to explain the ex- 
perimental elastic scattering cross sections. The GPDs 
can be studied in the ERBL region by such exclusive re- 
actions. 

Using the meson-pole model for the GPDs, we esti- 
mated the magnitude of the cross section of several re- 
actions pp — > pirN(A). We found a cross sections are in 
the range measurable at the high-energy hadron facilities. 
Experimental studies will be possible in two kinematics: 
one corresponding to the baryon B been slow in the rest 
frame, and another when B originates from the projectile 
(cross sections for the two kinematics are equal and we 
presented only one in our plots) . Experimental investiga- 
tions in the first kinematics would require both forward 
and recoil detectors, while in the second kinematics a 
forward multiparticle detector may be sufficient. Our re- 
sults indicate that both the nucleonic and the V — * A 
transition GPDs could be measured in hadron reactions. 
Future high-energy hadron facilities such as J-PARC and 
GSI-FAIR will provide opportunities to investigate var- 
ious interesting aspects of the GPDs in a way comple- 
mentary to the measurements using lepton scattering. 
These measurements will be valuable for finding the ori- 
gin of the nucleon spin, especially on the orbital-angular- 
momentum contribution. Note that the studied GPDs 
probe the isovector part of the quark angular momenta 
both in the p —> n and N — > A channels. 
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